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ABSTRACT 

We investigate the survival of planetesimal discs over Gyr timescales, using a unified ap- 
proach that is applicable to all Keplerian discs of solid bodies - dust grains, asteroids, planets, 
etc. Planetesimal discs can be characterized locally by four parameters: surface density, semi- 
major axis, planetesimal size and planetesimal radial velocity dispersion. Any planetesimal 
disc must have survived all dynamical processes, including gravitational instability, dynam- 
ical chaos, gravitational scattering, physical collisions, and radiation forces, that would lead 
to significant evolution over its Ufetime. These processes lead to a rich set of constraints 
that strongly restrict the possible properties of long-lived discs. Within this framework, we 
also discuss the detection of planetesimal discs using radial velocity measurements, transits, 
microlensing, and the infrared emission from the planetesimals themselves or from dust gen- 
erated by planetesimal collisions. 

Key words: planets and satellites: formation - Kuiper Belt - minor planets, asteroids - grav- 
itational lensing - Solar system: formation - stars: formation 



1 INTRODUCTION 



Terrestrial planets and the cores of giant pl anets are generally bel ieved to have formed hierarchically: small solid bodies ('planetesimals') 
condense from the gaseous ci rcumstella r disc I Johansen et al .1 120070 . collide repeatedly, and accumulate into larger and larger assemblies 
I Safronov|[l97i: lOoldreich et al. 2004a; iReipurth et al. 200^. Several aspects of this complex process remain obscure, in particular (but 
not limited to) the formation of planetesimals fr om dust grains i Blum & Wurm 200 sh . how to grow Uranus and Neptune in the short time 



available before th e gaseous disc is dissipated d Ooldreich et al.l 120 04 Jbh? t he origins of the large eccentricities of the extrasolar planets 
I Tremaine & Zakam ska 2004), the role of planetary migration ( Goldreich & Tremainjl 19801 ; Papaloizou & TerquemlboO^ . how the residual 
planetesimals wer e cleaned out of t he solar system i Goldreich et al.l2004bh . and why the solar system is so different from known extrasolar 
planetary systems ( Beer et alj2004 ). 



Given these large gaps in our understanding, it is worthwhile to investigate not just the difficult question of how planets form but also the 
simpler question of whether they can survive once formed. Planetary systems in the Galaxy have presumably been formed at a more-or-less 
constant rate, so it is reasonable to assume that most of today's planetary systems are at least several Gyr old. They must therefore have 
survived all dynamical processes - gravitational instabilities, collisions, viscous stirring or two-body relaxation, etc. - that would lead to 
a substantial change in their properties on timescales less than about 3 Gyr. Understanding what long-lived planetary systems are possible 
should help us to understand to what extent the properties of actual planets are shaped by the formation process as opposed to evolution 
('nature versus nurture') and may guide observers in searching for novel types of planetary systems. 

Orbiting solid bodies are often called 'planetesimals', 'planetary embryos', or 'planets' depending on their mass, but for simplicity we 
shall use the term 'planetesimal' to describe any body, whether solid like a terrestrial planet or gas-dominated like a giant planet, that is large 
enough so that gas drag and radiation effects (Poynting-Robertson drag, radiation pressure, Yarkovsky effect, etc.) are negligible. 

We shall find it useful to classify discs as 'hot' or 'cold' depending on whether or not the planetesimal orbits cross. Within the solar 
system, the planets form a cold system (except for Pluto) while the asteroid and Kuiper belts are hot. Among hot planetesimal discs, an 
important special case is the 'collision-limited' disc, in which the collision time between planetesimals is equal to the age of the disc. 
Collision-limited discs are likely to arise from discs in which there is a distribution of planetesimal sizes: smaller planetesimals have shorter 
collision times and therefore are destroyed first, so the dominant planetesimal population (by mass) always has a collision time that is roughly 
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Table 1. Sample planetesimal discs 
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Note: m.™/,^ (equation 1651 ). JTiJ^fi^ (equation 1741 ) are the minimum planetesimal mass for cold and hot discs, respectively; 
"^max" maximum mass for warm discs (equation 1751 ). A^^^J^ (equation 1671 ) and N^^-^ (equation 1741 ) are the 

maximum number of planetesimals per octave for cold and hot discs, respectively; N^f^^ (equation 1751 ) is the minimum 
number for warm discs. 

equal to the disc age. We shall also use the term 'warm' to describe discs in which the planetesimal orbits cross, but the impact velocities are 
so low that the cumulative effect of collisions does not substantially damage the planetesimals ( ^3.2.21 . 

The most important observational signature of many planetesimal discs arises from dust formed in recent planetesimal collisions. These 
'debris discs' were first detected from the thermal emission of the dust, which creates an infrared (IR) excess in the spectral energy distribution 
of the otherwise normal stars that they surround l Aumann et alj[l984 ). Spatially resolved debris discs are sometimes also visible from their 



scattered light. The IR excess (bolometric) luminosity relative to the lumi nosity from the pa ren t sta r depends on the distance and spectral 
type of the host star as well as the disc radius, but is typically > 10~^ (see lZuckermMlboOll and lwvatfcoOSi for reviews). The asteroid and 



Kuiper belts in our own solar system can be thought of as debris discs, although they would not be detectable around other stars with current 
technology since the bolometric IR excess is ~ 10"'^ for both belts. 

Debris discs have been detected around stars with a wide range of spectra l types (A to M) and ages (~ 10^ to 10^'^ yr|^ The dust 
masses inferred from these observations are 10"^ < Afdust/Me < 1 (Fig. 3 of lwyattlboosh . although this result depends on the assumed 



size distribution of the dust. The term 'debris' emphasizes that the lifetime of the dust grains from grain-grain collisions, Poynting-Robertson 
drag, or radiation pressure is robustly and considerably less than the stellar age, so the dust cannot be primordial and must be continuously 
regenerated, presumably by ongoing planetesimal collisions. 

We stress that the discs we consider in this paper are far more general than debris discs: they include hot discs in which the rate of dust 
generation may be undetectably small, cold discs in which there are no collisions, planetary systems, asteroid belts, planetary rings, etc. 

We begin by constructing a simple model for a planetesimal disc in [j2] In ij3] we describe the dynamical processes that act on planetes- 
imal discs. Collision-limited discs, which are a special subset of hot discs, are described in [j4] Non-gravitational forces on dust are briefly 
reviewed in [J5] The properties of long-lived discs are discussed in [jS] where we also study six sample discs (Table [Hi. We discuss possible 
techniques for detecting and studying planetesimal discs in [jT] and we summarize and discuss our results in [Js] 



2 A SIMPLE MODEL FOR A PLANETESIMAL DISC 

We consider a system of planetesimals orbiting a solar-type star of mass Mq . The extension to other types of stars is straightforward, 
but at this stage adds excessive complication. We shall focus on discs with an age to — 3 Gyr since these are likely to be more common than 
y ounger discs. We assume that the discs are gas-free, which is consistent with observations for most discs older than a few Myr (see Fig. 2 
of lWvatfcoOSi) . 



The surface density of the disc at semi-major axis a is written E(a); more precisely, the mass with semi-major axes in the range 
[a, a + da] is dA/disc = 27rEa da. The orbital period is 27r/r2, where 



In most cases, we shall assume that this material is collected in identical spherical bodies of density pp, radius r, and mass m — 
4:-Kppr'^/3 - a monodisperse planetesimal system (see the end of i]4.3l for further discussion of this approximation). The densities of the 

1 IR excesses characteiistic of debris discs have also been detected around white dwarfs, and they appear to be strongly coirelated with metal contamination 
in the white-dwarf photosphere, presumably arising from accreted planetesimals ( .Farihi et al..2009.) . 
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planets in the solar system range from 5.5 g cm~^ (Earth) to 0.7 g cm~^ (Saturn); in general gas-giant planets have smaller densities than 
terrestrial planets, but we shall sacrifice accuracy for simplicity and assume that all planetesimals have pp = 3 g cm"''. The surface number 
density of planetesimals is 

N=-. (2) 

m 

We assume that the eccentricities e and inclinations i of the planetesimals follow a Rayleigh distribution, 

cFn (X eiexp [— (e/eo)^ — (i/io)^] dide, (3) 

where eo and io are the root mean square (rms) eccentricity and inclination. The density of the disc in the direction normal to its symmetry 
plane (which we call the z-direction) is given by 

dM = n{z) dz = noexp{—^z^ /h^) dz, (4) 
where h = aio /-^2 is the rms ^-coordinate. The midplane number density no is related to the surface number density by 

m = —f=T = —j= — -■ (5) 

The radial velocity dispersion is 

f2aeo 



- ^2- 

The ratio io / eo can in principle have a wide range of values depending on the dynamical history of the disc, but in a variety of theoretical 
models and observed astrophysical discs io/eo ~ 0.5 jPones & Tremainjl993l : lstewart & Idj2000l) so we shall adopt this value throughout 
the paper (see also discussion following equation 1291 ). 

With these assumptions, the local properties of the planetesimal disc are specified by four parameters: the semi-major axis a, the surface 
density E, the planetesimal radius r or mass m, and the rms eccentricity eo. (We take the planetesimal density pp, the disc age to, and the 
stellar mass Mq to be fixed.) Collision-limited discs are specified by three parameters, since the requirement that the collision time equals 
the age provides one constraint on the four parameters. 

We shall find it useful to express many of our results in terms of the following three dimensionless parameters: 



Mo' 

ppS^ 

Mc7 

,2 



-=-»-'«'(3i^)(ioiur' 



The parameter p is a dimensionless mass, of order the ratio of the disc mass to the stellar mass. The parameter is a dimensionless density, 
which is almost always quite large for planetesimal discs, so we may assume u ^ 1 when necessary. The parameter r is a dimensionless 
age, scaled by Q and by Ea^/A/o since most evolutionary processes run slower if the orbital time is longer or the disc contains less mass. 

We parametrize discs in terms of surface density E and semi-major axis a since these are likely to vary much less than planetesimal 
mass m or rms eccentricity eo during the planet formation process. As points of reference, we shall refer to a set of six discs (Table[T} with 
semi-major axes a — 1 AU, a = 10 AU and a = 100 AU (a typical radius for spatially resolved debris discs), and with dimensionless mass 
p = lO"** and 10~^. The former mass corresponds to Mdiac ^ ttEo^ = TrpM p ~ O.SMjupitor — lOOM®, about the solid mass needed to 
form the giant planets and comets in the solar system ( Goldreich et al J2004b ). while the latter mass corresponds to Mdisc — M^. 



3 DYNAMICAL PROCESSES IN THE PLANETESIMAL DISC 

There are a few general criteria that planetesimal discs must satisfy. Firstly, we require that the planetesimal eccentricities and inclina- 
tions are not too large, which can be interpreted as the condition that most of the planetesimals are bound (e < 1), or that the disc is thin 
(h < a), or that the radial velocity dispersion is less than the circular speed (ar < fla). We shall write this criterion as 

eo<fe or -<0.35/e or ^ < 4t- W 
a ila v2 

In this paper, we adopt = 0.5. (See Table|2]for a summary of the dimensionless numbers used in our study.) 

Naturally, non-trivial planetesimal discs must contain more than one body. Assuming that the surface density does not vary strongly 
with radius and the disc is not too extended, the disc mass may be written 

Mdisc = 7r/„Ea^ (9) 

with fm of order unity, so the criterion that the number of planetesimals A*' — Aldisc/m > 1 becomes 

m < nfr^Ea^ or nf^Afa^ > 1. (10) 
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Table 2. Summary of dimensionless numbers 



Symbol 


Assumed value 


Defining equation and/or section 
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03, AppendixE] 
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(27]. Appendix|A] 
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22.67 


(|30), Appendix|A] 
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12.94 


go), Appendix|A] 


k 


1 
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0.46 


(56) 


M 


various 


m 


V 


various 


m 


T 


various 


m 


a 


-0.13 




b 


0.44 


ED 



In this paper we adopt fm = 1, which corresponds to a disc with a radial width of about 0.5a. 
We shall also require that the total disc mass is small compared to the stellar mass, 

where we adopt fd = 0.3. 

The total cross-sectional area of the disc is nr^N; thus, in the absence of mutual shadowing the geometrical optical depth of the disc as 
seen from the host star is 

r.-^- ^-^M.r\ (12) 
Note that this is differs by a factor of order unity from the normal geometrical optical depth of the ring, A/'vrr^. 



3.1 Cold discs 

A 'cold' disc is one in which planetesimal orbits do not cross. A planetesimal in an orbit with eccentricity e has a total radial excursion 
of 2ae. The typical radial separation between planetesimals is 

Thus most orbits do not cross if the rms eccentricity is 

-o<-^>^ot = /c^, (14) 

and we shall choose /c = 1, at which point just over half of the particles cross if their semi-major axes have a Poisson distribution. (A 
refinement of the preceding condition is to include the radius of the planetesimal in the crossing condition, that is, 2(ae + r) < Aa, but in 
the cases of interest to us this correction is unimportant.) Equivalently, one can describe cold discs as those for which Gt < Uhot, where 

_ fc mn _ fc ^/UM^py 



3.1.1 Gravitational stability 

Gravitational instability in cold discs is associated with dynamical chaos, which leads to growth in the eccentricities and inclinations of 
the planetesim als. The growth rate can be extremely slow - for example, tens of Gyr for Mercury in the current solar system configuration 
i Laskajliooi) . No rigorous analytical formulae for the rate of chaotic evolution in cold discs are available. However, N-body experiments 
suggest that cold discs can survive for millions of orbits if the separation l ll3t is typically a few times larger than the Hill radius 
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A / \ 1/3 



Thus, stability requires 

^> Halt J " Mo>'-'^'"^"^- 

Chambers et al. ( 19961) find k — 7-1 1 from integrations of planetesimals with masses between 10~^ and Mq lasting ~ 10^ yr[f| They 



also find that the factor k varies slowly with planetesimal mass, so a stability criterion that fits their results more accurately is 

a \ Mq J Mq 

This result depends on the duration of the integrations, but only weakly: increasing the duration by a factor of ten typically increases the 
minimum stable separation by about one Hill radius. 

The criterion ilSi was derived from simulations with planetesimals of mass 1O~^-1O~^M0; the extrapolation to larger planets is 
somewhat uncertain but probably not a major source o f error. For J upiter-mass planets the formula predicts Aa/a > 0.5, corresponding to 



fc ~ 7 in the notation of equation l ll7t . For comparison. I Juric & Tre maine (2008) estimate k ~ 12-14 from orbit integrations of planets with 



masses between 0.1 and 10 Jupiter masses; their larger value of k probably arises because they used a range of masses rather than a single 
common mass for the planets. In any event, equation J18t should be correct to within a factor of two or so. 



When the planetesimals in a cold disc have non-zero eccentricities, larger separations are required for stability (e.g.,|Yoshinaga_etal 



Il999h : however, we shall not include this refinement since the detectability of cold discs does not depend on the rms eccentricity (equation 

[Ml). 



3.2 Hot discs 

A hot disc is one in which planetesimal orbits cross. In hot discs, the rms eccentricity or radial velocity dispersion exceeds ehot or Uhot 
respectively (equations 1141 and 1151 ). As we discuss below (equation 1401 ). 'warm' discs may also be defined, in which the orbits cross but 
the collisions do not substantially damage the planetesimals over the lifetime of the disc. 



3.2.1 Gravitational stability 

The disc must be gravitationally stable to axisymmetric perturbations i Toomre|[l964 ). The Toomre stability criterio n is derived from the 



WKB dispersion relation for axisymmetric density waves. For a fluid disc in a Keplerian potential this may be written as iBinnev & Tremaine 
2008h 



where uj is the frequency, A is the radial wavelength, and 



-GE V2.,' ^20) 
Acrit = — ^ — = 47r a/i. 

This dispersion relation was derived for a barotropic fluid whereas the planetesimal disc more closely approximates a collisionless fluid, but 
for our purposes the results should be accurate enough. 

The disc is stable to perturbations at a given wavelength A if cj^ > 0, which requires 

Q^>4-^(l--^). (21) 



Acrit \ Acrit , 

The maximum of the right side occurs at A = Acrit /2 and equals unity, so the disc is stable to perturbations of all wavelengths if Q > 1. 
However, this result needs to be modified if the number of planetesimals is so small that the disc cannot be approximated as a continuous 
fluid. The typical separation between planetesimals is given by equation l ll3t . so the number of planetesimals in one wavelength A is N\ = 
A/Aa = 27raEA/m. The continuum approximation should be valid if N\ ^ 1 or A > Ac = fqin/ {2-naT,) with /q of order unity. We 
adopt /q = 1 when a choice is necessary. Then a necessary requirement for stability is that equation is satisfied for all A > Ac, or 



q2 I - " "c ^ -.cm/ ,22N 

^ ^ 4(Ac/Acrit)-4(Ac/Acrit)' if Ac > Acrit/2. 



1 if Ac < Acrit/2, 



Note that lChambers et alj|l996l) define the Hill radius as (2m/3M0 )i''^, a factor 2^/^ larger than our definition. 
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cold 




m/M3 



Figure 1. Stability and other properties of planetesimal discs, as a function of planetesimal mass m and rms eccentricity eg. Tlie parameter ^ = T,a^ /Mq ~ 
A/disc/Af© is approximately the disc mass relative to the stellar mass. The region to the right of the green, vertical solid line, labelled 'N < 1', is not 
allowed because the planetesimal mass exceeds the assumed disc mass. The slanted solid line divides hot discs, in which planetesimal orbits cross, from cold 
discs (equation 1141 ). The regions labelled 'unstable' represent hot discs that are unstable (in pink) and cold discs that are unstable (in blue) according to the 
approximate equations derived at the end of ^3.2.11 The slanted dashed line separates hot discs in which the encounters are shear-dominated from those in 
which encounters are dispersion-dominated. The pink and blue shaded regions represent the allowed parameter values for hot and cold discs, respectively. 



This may be rewritten as 



eo 



> V2-KI1. if — — < 



Mp 



fQ 



eo > 



1 - 



Otherwise; 



(23) 



there is no constraint if the square root in the second equation is negative. 

The relation between the stability criteria for hot and cold discs can be clarified using Figure[T] which plots allowed regions as a function 
of the dimensionless parameters m/M© (ratio of planetesimal mass to stellar mass; horizontal axis) and eo (rms eccentricity; vertical axis). 
The axes are logarithmic. We aim for clarity at the sake of accuracy by neglecting all factors of order unity for the rest of this subsection. The 
diagram shows the following constraints: 

(i) There must be at least one planetesimal in the disc (equation 1101 ). so m/Ea^ < 1 or m/A/© < /i where /i is defined in equation 
(|7}. The boundary m/M© = ^ is represented by a green, vertical solid line; the excluded region to the right of this line is labelled 'N < V . 

(ii) The division between hot and cold discs ( 114b may be written eo < m/ ^Mq, which is marked by a slanted solid line. 

(iii) The condition (TsJ for the gravitational stability of cold discs is m/MQ > fi^'"''^, which appears in the figure as a vertical dotted 
line; the unstable region to the left of this line is labelled 'unstable' (in blue). 

(iv) The condition J23t for gravitational stability of hot discs becomes eo > A* 'f m/MQ < . The second of equations ( I23t is 
neglected because it applies only over a range of a factor of two in planetesimal mass. The unstable region is bounded by a solid horizontal 
line and labelled 'unstable' in pink. Although this derivation was carried out for hot discs, it should apply to cold discs as well so long as 



m/M(i 



however, it adds no new restrictions on cold discs since since these are already unstable by condition U8b . 



The pink and blue shaded regions represent the allowed parameter values for hot and cold discs, respectively. 
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3.2.2 Collisions 

The collision time in a system with isotropic velocity dispersion a and number density n is given by i Binnev & Tremainel2008 , equation 



[7.195]) 

= 160Fn(jr^ (1 + O) , (24) 

where 

e - ^ (25) 

is the Safronov number. The factor (1 + 0) reflects the enhancement in the collision rate due to gravitational focusing. If the velocity- 
dispersion tensor is anisotropic, we define 6 by replacing cr by the radial dispersion cr,., and then with equation ^ we have 

, / X 2/3 1/3 



Merel \MqJ eg 

where the dimensionless density is defined in equation 

Equation l l24b requires several corrections for use in discs. Firstly, the velocity-dispersion tensor in a Keplerian disc is not isotropic; for 
the Rayleigh distribution l|3j and our choice io / eo — 0.5, this requires replacing a with the radial dispersion ar and the factor (1-1-0) with 
(/i -|- /2O) where /i = 0.690 and /2 — 1.521. Secondly, we must account for the falloff in density above the disc midplane, which we 
do by replacing the number density n in equation i24i with the number-weighted average J n^{z)dz/ J n{z)dz = 7io/\/2, where n{z) is 
given by equation (O. Thus, we replace equation l l24t by 

t-' = 2''^\^^\oay (/i + /2e) 

(27) 



for_io/eo ~ 0.5, as we assume throughout. See Appendix |A] for details on how to modify the more general formulae of IDones & Tremaine 



1 19931) to arrive at the preceding result. 

Equation ( 127b neglects the Keplerian shear in the disc and thus is only valid for dispersion-dominated encounters, for which 

(Tr > max {r, rn} , (28) 

where the Hill radius ru is defined in equation ( I16t . Note that ru > r if the dimensionless parameter u (equation (T)) exceeds 0.7, which is 
almost always true, so in practice condition ( I28t reduces to 

cTr > firn or eo > (m/M©)^''^ (29) 

This boundary is marked as a slanted dashed line in Figure [T] The allowed region in parameter space below this line, in which encounters 
are shear-dominated, is relatively small but still requires consideration. Note that condition J29t can be rewritten with the help of i26l as 
e < so shear-dominated encounters occur only if the rms eccentricity is so small that the Safronov number is greater than unity by 
the large factor u^^'^. When the condition in equation l |29b is violated, the disc is so flat that encounters may excite out-of-plane motions less 
efficiently than in-plane motions, so that io /eo may be less than our assumed value of 0.5. We do not, however, attempt to model variations 

in io / eo in this paper. 

Formulae for the shear-dominated collision rate are given bv lOreenzweig & Lissauei 1992 ) and Pones & Tremaine 1993h . Adapting 



these formulae to the present model and notation (see Appendix IaI. we have 

-1 _ f^AAfi^aVa-i {m/MQf^-^ , ^tyiIv^I^ < < f^rn, 

' [mnr'/^a'/^m/MQy/\ ar<nrn/^^'\ ^^^^ 

where fi = 22.67, /s = 12.94, and v is defined in equation l|7j. 

Collisions between equal-mass planetesimals may have various outcomes. If the gravitational binding energy of the planetesimals is 
negligible compared to their relative kinetic energy (Safronov number B ^ 1) then (i) high-speed collisions will shatter the planetesimals 
and disperse the fragments, while (ii) low-speed collisions will leave the planetesimals unaffected or produce small craters. If the gravitational 
binding energy is much larger than the kinetic energy then (iii) high-speed collisions will still shatter and disperse the planetesimals, while 
(iv) the outcome of low-speed collisions will be a single gravitationally bound object containing most of the mass of the two planetesimals 
(see ij4T|for further detail). For our purposes, any of the outcomes (i), (iii), or (iv) leads to a substantial change in the mass distribution of 
planetesimals if the collision time is less than the age of the disc, tc < to - and therefore is inconsistent with our requirement that the disc is 
long-lived in its present state. 

However, many collisions of type (ii) could occur without substantially altering the planetesimal mass distribution. We therefore differ- 
entiate between dynamically hot discs, in which the collision time is longer than the age, tc > to, and 'warm' discs, in which the collision 
time is shorter than the age but the velocity dispersion ar is small enough that the collisions do not substantially damage the planetesimals 
over the lifetime of the disc. The survival criteria for warm discs are discussed in ^3.3\ 
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3.2.3 Gravitational scattering 

The timescale for a substantial change in the rms eccentricity and inclination due to gravitational scattering of planetesimals must also 
be longer than the disc age. We write this as tg > to, where 



dlog eg 
dt 
fim 



[Si(io/eo) C + S2(io, eo,m/MQ)] . (31) 

en 



Here the term S'i(io/eo)C incorporates the effects of close or dispersion-dominated encounters, in which the relative velocity is dom- 
inated by the velocity dispersion cr,., while 6*2(10, eo, tti/Mq) represents the effects of distant or s hear-dominated encou nters. The factor C 
is usually called the Coulomb logarithm. We evaluate these terms by specializing the formulae of Stewart & Idal 1 2000l) to the case where 



the colliding planetesimals have the same mass and the same eccentricity and inclination distributions. From equations (3.29), (3.31) and 
(6.6)-(6.8) of that paper, we find (0.5) = 4.50 and 

2C = log(A' + 1) - log(A^ + 1) + 



A2 + 1 A? + r 



M© 2r 2 , .2x 

Ac = (eo + io) 

m a 



1 + 



M0r(e2 + i2 + i(|m/Me)2/3) 
From equations (4.8) and (6.5) of lstewart &"ldal ilOOd) . 



1/2 

(32) 



52(io,eo, m/Me) = ^'f ^ [W{e/el) - W{e/il)], (33) 

where e= l^{lm/Mef^\ 



exp(r-y) 



X (34) 

y 



and ^ ~ 2 is determined by an empirical fit to N-body simulations. 
Using equation l l26t , equation | |3U may be rewritten as 



= AAf^r-^e^SiC', (35) 

where 

C' =C + S2/S1 (36) 

is the correction factor arising from the Coulomb logarithm and shear-dominated encounters. 

The characteristic timescales l |27| l and l |35b for physical collisions and gravitational scattering can be combined into a single relaxation 
timescale, 

i7cL = ic'' + = mmr^ (/i + /2e + /sC'e") , (37) 

where fs — Si/16 = 0.28. A long-lived disc must have trciax to. 

The effects of collisions and gravitational scattering can be clarified using Figure[2l which plots allowed regions as a function of m/M© 
and eo, in logarithmic coordinates. We neglect all factors of order unity for the rest of this subsection. In this approximation, the requirement 
for survival of hot discs simplifies to 

This result holds only if the collisions are dispersion-dominated, but if they are shear-dominated the relaxation time becomes shorter than 
this formula would predict so the criterion l |38t remains necessary (but not sufficient) for survival of the disc. 

The results can be re-written in terms of the dimensionless time r defined in equation l[7}. Discs with Safronov number O > 1 have 
m/MQ > Col \/v and the survival criterion OSt is m/MQ > jv^ (6 < 1) or m/MQ < eg/r {Q > 1). These constraints are shown 
in Figure[2] The excluded regions are bounded by the cyan, vertical solid (m/MQ — r'^/u'^) and green, slanted dot-dash (eg = mr/Mo) 
lines, and are labelled 'trciax < to' - Notice that for a given value of the rms eccentricity eo or velocity dispersion ar, there is only a finite 
range of planetesimal masses or radii in which the relaxation time exceeds the age: the disc cannot survive if the planetesimal mass is either 
too small or too large. This behavior arises because for small masses m or radii r (0 <C 1), we have froiax oc r oc m^^^ at a given surface 
density and eccentricity, while for large masses (0 3> 1), we have treiax oc oc r~^. 
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Figure 2. Survival of hot discs as a function of planetesimal mass m and rms eccentricity eg. The parameters t and u are defined in equation (7). The slanted 
solid line divides discs according to their Safronov number ©. The regions marked by 'trelax < to' are not allowed because the collision time (when © < 1) 
or the gravitational scattering time (when O > 1) is shorter than the disc age. The pink shaded region represents the allowed parameter values for hot discs. 



3.3 Warm discs 

As described in i ]3.2.2l planetesimals in warm discs may suffer collisions but the velocity dispersion is low enough that these collisions 
do not destroy the planetesimals. Thus a necessary condition for the survival of warm discs is that the squared velocity dispersion a'^ must 
be less than Q^, the energy per unit mass required to disperse the planetesimal into fragments, which we obtain from equation ( 144b below. 

This criterion is not sufficient, for two reasons. Firstly, if the gravitational binding energy of the planetesimals is much larger than 
the kinetic energy, an inelastic collision is likely to leave the colliding planetesimals as a gravitationally bound pair, which alters the mass 
distribution and therefore is inconsistent with our assumption that the disc has not evolved. Therefore we require that the Safronov number 
O < 1 for warm discs. 

Secondly, low-velocity collisions can chip or crater the planetesimals even if they are not disrupted. Thus if the collision time tc is 
much less than the age to the planetesimals may be gradually eroded away even if they are not dispersed in a single collision. A simple 
para metrization of the erosive pro cess is to assume that the mass lost in a typical collision of two objects of mass m at relative velocity v is 
(e.g., lThebault & Augereai]|2007l) 



Am = 0.5m (^—j , < Q*d (39) 

(the factor 0.5 arises because the usual definition of a dispersive impact is one in which the mass of the largest fragment is less than half of 
the target mass; see ^4A\ . 

Erosion through impact involves a number of complicated processes such as crack propagation in brittle materials and plastic flow 
in ductile materials, and melting or sublimation at high impact velocities. Our present state of knowledge is derived from a variety of 
approximate physical models, experiments, and numerical simulations, most intended for situations far removed from planetesimal collisions. 
These typically yield values of 7 between 1 and 1.5. For example, (i) impacts of small pellets into targets composed of ice-silicate mixtures at 
speeds up to 12 km s~'^ yield 7 ~ 1.2 IjC oschny & Griin 2001); (ii) numerical simulations of collisions between rocky and icy bodies yield 



yield ^ 

7 ~ 1 1 Benz & Asphaug 19991 : Stewart & Leinhardtl2009 ). We shall adopt 7=1 recognizing that this is (a) oversimplified; (b) conservative 



(in that the actual erosion rate is expected to be smaller for low-velocity collisions if 7 is larger). 

The criterion for survival is that the cumulative mass loss Am{to/tc) < m. Thus warm discs must satisfy 

tc <to, e < 1, and < (tc/to) Qd- (40) 
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where is given by equation iAAV\ 

An additional constraint for warm discs is that collisions do not result in excessive 'viscous' spreading of the disc. A disc initially 
localized at semi-major axis a that loses a small amount of energy 5E must spread by 5a <C a (,Brahic,1977,) . where 

If the typical energy lost per unit mass in a collision is /eo-^, where /e ~ 1, then the rate of energy loss is faa'^lMdiac/tc- Therefore, in a 
time to the disc spreads to 



Sa / 32 fQaa^ to 



a y GMq tc 

Requiring Sa/a < 1 and assuming B <c; 1 yields 



(42) 



fe J \0.69 J V3gcm-M05cmy VlO-^'SGyr/ VlOAU/ 
This constraint does not restrict the allowed range for warm discs among the six sample discs considered in this paper, i.e., it is satisfied so 
long as the warm discs satisfy all of the other constraints we have already discussed. 



4 COLLISION-LIMITED DISCS 
4.1 The effects of collisions 



Collisions or impacts may crater, shatter or disperse the target, as discussed in i ]3.2.2l and H3.3\ The distinction betw een shattering and 
dispe rsive impacts arises because the pieces of the target may remain gravitationally bound even after the target is shattered ISenz & Asphau j 
1999h . 



The outcome of an impact depends mainly on the relative speed, masses, and composition of the impacting bodies. Consider an impact 
between two planetesimals of masses mi (the target) and m2 ^ mi (the projectile), at relative velocity Av. The kinetic energy of relative 
motion is Ek = ^rfiAv^ where m = m\m2/{m\ + m2) is the reduced mass. The impact is shattering if the kinetic energy per unit target 
mass Ek/m\ > QJ, and dispersive if Ek/mi > QJ,, where Q* is a function of the mass and composition of the target. Roughly speaking, 
Q* (m) is the binding energy per unit mass of the target - the energy needed to rupture its internal chemical bonds for Qg and this energy plus 
its gravitational potential energy for Q^. In the strength-dominated regime, when the target is small enough that its self-gravity is negligible, 
we expect Qg ~ QJj, and both should be independent of m in the idealized case in which the strength of the target is independent of its size. 
In practice, both simulations and experiments find that Q* declines slowly with m in the strength-dominated regime. In the gravity-dominated 
regime, Qg <^ Q*q (an impact that shatters a target may not impart escape speed to the fragments, so they reaccumulate as a gravitationally 
bound rubble pile) and we might expect that Q oc m?^^ since the gravitational binding energy per unit mass of a homogeneous body scales 
as m?^^. 

The dependence of Q*jj on target mass is typically parametrized as 

ODM = Qo(a;S + 4), (44) 
where xo = m/mo; Qo, mo, a and h are parameters to be fitted to experiments or simulations. We shall adopt mo = lO^'^ g, Qp 



6 X 10 ^ erg g ^, a = —0.13 and b = 0.44, these values being compromises between the results for ice and basalt given by Benz & Asphaud 



I I999I) . This simple form fails for very small mass, since it pr edicts Q*n{rn) — >■ 00 as m — > 0; experiments with high-velocity impacts of 



small bodies suggest QJj « lO'^ erg g ^ independent of mass jplynn & Durd j2004l) . so for a;o ^ 1 we use the smaller of this value and the 
prediction of equation l l44t - the transition occurs at m ~ 10** g. 



4.2 The coUisional cascade 



When the collision time is much less than the age of the disc, the mass distribution of the smaller planetesimals is likely to be established 
by a 'collisional cascade', in which large bodies are dispersed by collisions into smaller bodies, these in turn being dispersed into smaller 
ones, until bodies of size < 1 /xm ('dust grains') are removed by Poynting-Robert son drag and radiation pressure (see jBll. We now derive 
an ap proximate form for the steady-state mass distribution in a collisional cascade jDohnanvi|[l969l : lo'Brien & Greenber J2OO3I : |Pan & Saril 
200a), using the following assumptions: 



(i) The number density of planetesimals per unit mass is a power-law function of mass, at least over a limited range, 

dn 

- — oc m . 
dm 



(45) 



The distinction between 'hot' and 'warm' discs based on whether or not the collisions are destructive is moot when the collision time is longer than the age. 
Our (arbitrary) convention is that such discs are 'hot', not 'warm'. 
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(ii) The rms eccentricity of the planetesimals is independent of mass. (This assumption is made for simphcity, but appears to hold approx- 
imately in the asteroid and Kuiper belts.) 

(iii) The velocity dispersion ^ is large enough that a typical impact between two bodies of equal mass has much more than enough energy 
to disperse the two bodies, i.e., al > Q^- which requires (Tr > 3 x 10^ cm s~^ (QJ,/10^ erg g"^)^''^ or eo » 0.005(a/10 AU)^''^ x 
(Q^/lO^ergg-^)^/^ 

(iv) The specific impact energy required for a dispersive impact is a power-law function of mass, Q J, oc . According to equation J44t . 
we expect j = a for m <^ mo (strength-dominated) and j = 6 for m 3> mo (gravity-dominated). 

(v) The cross-section for a collision between two bodies of masses mi and m2 ^ mi is (Tcoii(mi, 1712) oc m'l. From equation i24i we 
expect / = I for 9 ^ 1 (ctcoii oc oc m^/^) and / = | for 6 S> 1 ((Tcoii oc r^Q oc rm oc m"*/^). 

According to assumption (iii), the smallest projectile mass that can disperse a target of mass m in an impact is mmin i'm) ~ mQ}^ ja'l <C 
m. The rate of collisions in which planetesimals of mass > m are dispersed is approximately 

^ , . , . dn(mi)dn(m2) , , 

<P„(m) oc / dmi / dm^ — ; ; ctcoii (mi, m2) . (46) 

This formula is accurate to a factor of order unity only, since it neglects the fact that one object of mass > m is dispersed when mi > m 
and m2 < m, while two are dispersed when mi , m2 > m; it also neglects the possibility that one or more of the collision fragments is more 
massive than m. These inaccuracies do not affect our final result. 

Since most of the fragments in the collisions that dominate this rate will have masses < m, the mass per unit volume in planetesimals 
larger than m decreases at a rate given approximately by 

^ , ^ ^ . ^ pi dn(mi)dn(m2) , , 

q?m, (m) ^ mq?„(m) oc m / dmi / dm2 — ■ ; ctcoII (mi, m2) . (47) 

With assumptions (i), (iv), and (v), we find 

$„(m) ocm='+^+'-^f2+^'. (48) 
In a steady state, the mass flux (m) must be independent of mass, so 

With the parameters used in this paper, almost all dispersive collisions have O <C 1 so we can set Z = | and obtain dO'Brien & Greenberg 

(50) 



2ooa) 

6 + 3j 

Matching the power-law behavior for m ^ mo and m ^ mo, we have i Lohne et"al]|2008t) 



dn(m) _ J (m/mo)-("+^'^/(^+^''\ m > mo, 



dm ^\ (m/mo)-("+^")/(«+3"', m < mo. 

These expressions are only valid if (i) the collision time at mass m is short compared to the age of the disc; (ii) the typical impact velocity is 
sufficient to disperse a body of mass m, cr,^ ^ Qni'T'-)' (iii) the velocity dispersion is independent of mass; (iv) the mass m is sufficiently 
large that the lifetime to radiation pressure and Poynting-Robertson drag is much larger than the collisional lifetime. 

For the exponents a — —0.13, b — 0.44 given after equation l l44t . we have p = 1.89 for m <^ mo and p = 1.68 for m 3> mo. 
With these exponents the total mass in the collisional cascade (oc J mdn cx m^~^) is dominated by the largest bodies, while the total 
cross-section (oc J m^^^ dn oc m'^''^^^) is dominated by the smallest bodies. 



4.3 Properties of collision-limited discs 

Collision-limited discs can arise if we assume that there is an initial distribution of planetesimal masses in which the total disc mass is 
dominated by small bodies, say dn oc dm with P > 20 In such discs the collision time is shorter for planetesimals of smaller mass; 
thus a collisional cascade is established for all masses below some maximum mmax. We have seen that the mass in the collisional cascade is 
dominated by the largest bodies in the cascade, while if P > 2 the mass in the 'primordial' regime m > mmax is dominated by the smallest 
bodies. Thus, the overall mass in the disc is dominated by bodies with mass ~ mmax, and for many purposes we may treat the disc as a 
monodisperse system composed of bodies with a single mass mmax. The difference from our earlier discussions is that now the planetesimal 
mass is not a free parameter; rather, it is determined by the condition that the lifetime of a planetesimal of mass mmax subject to dispersive 
impacts is equal to the disc age. We now derive this condition. 

* The condition P > 2 appears to hold for most planetesi mal populations in the solar system. For the classical and excited Kuiper belts, P ~ 3.3 and 
2.1 respectively, and for t rans-Neptunian objects P ~ 2.5 ifiernstein et alj|2004). Numerical m odels of the formation of the Kuiper belt give P = 2.3 
iKenvon & Bro mlev'2004). Asteroid observations yield P between 1.9 and 2.3 iParker et alj2008l) and Jupiter-family comets have P = 1.9 though with large 
uncertainties fpernandez 20 05) . 
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Following the discussion in the preceding section, we shall assume that dispersive impacts have Safronov number <C 1 and are 
dominated by impactors of mass much less than the target mass. Then if we equate the collision time from equation l l27b to the disc age to, 
we have 



f 3mn 



2/3 



dno (m') , 
. . am , 



dm' 



(52) 



in the midplane of the disc is 



we have 



/a^ is the smallest impactor that will disperse a target of mass mmax- Assuming that the mass distribution 



dno {m) rin 



dm 



tn 



5.32 /i 



(m/mn 
(m/m„ 



2/3 

^maxm-ixiax^ r 



m 5C m„ 

m > TTln 



2/3 



1-p 



p-1 0.69 

The total midplane mass density po is related to the surface density by S = ^/iipoaio (equation (5)) and po is given by 



(53) 



(54) 



dm + m„ 



^2-p 

assuming p < 2 < P. If we divide the second equation by the first and use equations ([5} and ^ we obtain 



nto 



1.25/7 



/i 
0.69 



1/3 2/3 



p-1 



(55) 



(56) 



which is an implicit equation for the characteristic planetesimal mass mmax. Here /? = (p — l)[l/(2 — p) + 1/(P — 2)]/5.32. We choose 
p = 1.68 and P = 4 so /r = 0.46; a steeper high-mass slope P = 6 would change this only to /t = 0.43. Alternatively, we may write 



1.25/7 



/i 

0.69 



-1 1/3 2/3 
m,niax Pp 



2(P-1) 



(57) 



By replacing the mass m in our discussion of monodisperse discs with mmax, most of the results of Sj2]and !j3]can be applied to collision- 
limited discs without major errors. For example: (i) The gravitational stability of hot discs depends on the surface density E oc J m dn, 
which is dominated by masses near mmax when p < 2 < P. (ii) For most purposes the appropriate replacement for the collision rate 
f^^ oc nr^(/i -I- /26) (equation J27I ) in a disc with a distribution of masses is the mass-weighted collision rate oc J mr^(/i + /26) dn, 
which in turn is proportional to J m^^^ dn for <C 1 and J mJ^'^ dn for G ^ 1. These integrals are dominated by masses near mmax when 
p < 8/3 < P and p < 10/3 < P respectively, (iii) The gravitational scattering rate is oc J r^Q^ dn cx J m^ dn so this is dominated 
by masses near mmax when p < 3 < P. All of these inequalities are satisfied for our nominal values p — 1.68 and P — 4. 

Finally we note an interesting inconsistency in the results we have derived so far. Neglecting factors of order unity, the collision time 
fc,m for a monodisperse disc with ^ 1 is given by equation ( 127b as t^Jn ~ noaj.m'^^^ / pp^'* . The analogous collision time tc,ci for a 
collision-limited disc is given by equation ([54} as t^^j ~ rimaxO"rmmax/Pp'''^(Qr)/''"r )^ ~ ic^m(Qr)/''"r)^ ?S ^c^m- The difference 
arises because the collision-limited disc has a large population of bodies with mmin ^ 'm ^ mynayi that can collide with and disperse the large 
planetesimals, whereas in the monodisperse disc these are destroyed only by collisions among themselves. What then is the state of a disc in 
which fc.ci ^ io ^ ic.m? Should it be regarded as monodisperse or collision-limited? We assume here that such discs are monodisperse but 
this assumption may be oversimplified. 



5 NON-GRAVITATIONAL FORCES ON DUST 

We also describe the most important non-gravitational forces on dust grains jBumsetal.ll 19791) . since the distribution of dust grains 
determines the infrared flux from debris discs; these forces can also be relevant for the planetesimals in warm discs. Gas drag on the dust is 
unimportant since we are focusing on discs older than a few Myr, at which point the gas in the protoplanet ary disc h as disappeared. The ratio 
of repulsive forces from the stellar wind and radiation pressure to the attractive gravitational force is (e.g.. lStrubbe & Chiana,2006i) 



Lr.,Vr 



Pp 



(58) 



IQ-K GMqcPpT - ■ ' \^3gcm-3y \^l/im 

in which we have assumed that the host star has the solar mass Mq and luminosity Lq, pp is the dust grain density, r is the grain radius, and 

Mvu 



— Qpr ' 



(59) 
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where M is the rate of mass loss in the stellar wind, V w is the wind speed, and Qpr is the radiation pressure efficiency factor (averaged over 
the stellar spectrum) as defined bv lBums et al. ||l97i. For stars with the Sun's luminosity and age, the contribution of the stellar wind to 



Vr is negligible (< 10 ^ ). Grains created b y collisions on circular orbits with the Keplerian speed are unbound if /3 ^ ^ - because their 
eccentricity is /3/(l — /3) jBurns et al.lll979l) - and thus P = \ defines the 'blow-out radius', 

rt = 0.38 /im Vr ■ (60) 



^3g cm 

The rate of orbital decay from Poynting-Robertson and stellar-wind drag is ( lBurnsetal.lll979l) 

1 _ Ida _ 3 LqV4, _ 2V4, GMq 



tpB. a dt Stt c^a^ppT Vr 

where 



P, (61) 



V4, = Qpr + (62) 
Lq 

Numerically, we have 

tpH = 8.0xl0^r/?-g(^)\ (63) 

The contribution of the stellar wind to V^ is negligible for particles with r > 1 pm but grows as for r < 0.3 /im, and equals the 
contribution due to Poynting-Robertson drag at r ~ 0.1 fim (Burns et al. 1979). 

The thermal emission from the dust is determined by its absorption efficiency factor Qa, which is similar in magnitude to Qpr and equal 
to it if scattering is neglected. In the geometric optics limit, Qa is independent of wavelength and close to unity for typical dust grains; in 
the limit of long wavelength, where X = 2nr/X ^ 1, Qa ~ X (equation 1921 ). The contribution of a grain to the thermal emissivity at 
wavelength A is proportional to ixr^Qa- Thus for a power-law mass distribution with exponent — p (equation 1451 ) in the range | < p < 2 (as 
is the case for a collisional cascade; see Sj4j2jf|the total thermal emissivity is dominated by grains with X = 2-Kr/\ ~ 1 or r ~ riR,(A) = 
A/(27r). More precisely, a fraction / of the emission comes from particles with radii that exceed rm./ k., where / = 1 — (3p — 5)«:^''''"^\ For 
p= 11/6, the expected value for a collisional cascade (we use p = 11/6 = 1.83 rather than p = 1.89 as at the end of i]4.2l for reasons given 
after equation 1941 ). 75% of the emission comes from particles with radii that exceed 0.25riR. For observations at A = 20 /xm, na = 3 pm 
and 75% of the emission comes from particles with radii that exceed 0.8 /xm, where /3 = 0.24 (assuming pp = 3 g cm~'' and Qa = 1, since 
most of the stellar emission is at shorter wavelengths). At longer observational wavelengths (3 is even smaller for the particles dominating 
the emission. Thus radiation pressure is negligible, except perhaps for accurate modeling at the shortest observational wavelengths. 

Poynting-Robertson drag is also negligible, at least for detectable debris discs, as shown by the following argument. Assume for 
simplicity that the dust particles have a single size. Using equations M2\ and J27| ( the collision time tc for dust is related to the geometrical 
optical depth as seen from the host star Tp by 

flt, = ^^ or = 3.6 X W yr frn ( ^] ^ (-^YU—^Y^^ ■ (64) 
64/irp ■' \0.69 J VlO-V VlOAU/ 

For particles with absorption efficiency Qa — 1, the o ptical depth Tp is equal to the bolometric luminosity of the disc relative to the star, 
which exceeds 10~^ in almost all observed debris discs ( Wvattl2008!) . Therefore the collision time 4 x 10* yr (a/10 AU)'^''^ in observed 
discs, shorter than the Poynting-Robertson drag time l |63l l, so the grains are destroyed by collisions before they experience significant orbital 
decay. 

Based on these arguments, we neglect non-gravitational forces on the dust distribution when calculating its emission properties. 



6 THE PROPERTIES OF LONG-LIVED DISCS 

The primary goal of this paper is to explore the properties of planetesimal discs that can survive for most of the age of the Galaxy. 
Even in our highly simplified model, the local properties of discs are specified by four parameters: semi-major axis a, surface density E, 
planetesimal radius r or mass m, and rms eccentricity eo, while collision-limited discs are specified by three parameters. It is challenging 
to visualize the properties of a four-dimensional parameter space. As a first step, we outline some general conclusions that arise from the 
discussion of the previous section. 



^ The assumption of a power-law mass distribution neglects the oscillations that appear in the mass distribution of a collisional cascade at radii that are not 
too far from the blow-out radius iKriv ov et aL20Q6.) ; these oscillations can change the number density at a given radius or mass by a factor of ^ 3. 
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6.1 Cold discs 

6.1.1 The minimum planetesimal mass in cold discs 

Eliminating the separation Aa between equations i\2>\ and dlSt yields the minimum mass in cold discs, 

m > m''^^ = IOOMq/x' *' ~ 60Me {-^) ' (65) 
Thus, the minimum planetesimal mass in a cold disc with = 10^* is 0.2 Jupiter masses, while if = 10"'^ the minimum mass is 0. 1 Earth 



6.1.2 The maximum number of planetesimals in cold discs 

If we write the disc mass as A/disc ~ T^fm^o? (with fm = 1 this is roughly the mass per octave in semi-major axis; a disc extending 
over multiple octaves in semi-major axis could have fm ^ 1), then with equation J65t we have 

iV = ^ < iV-':^ = % ('^'l , (66) 



^ ^ ■ max — o \ 

m 8 \ m 

or TV ~ 1 for Jupiter-mass planets, 6 for Earth-mass planets, and 20 for lunar-mass planets. 
This result can be re-cast in terms of the surface density, 

/ /; N -0.43 

N < iV;^;^ = (y^) • (67) 

Thus cold discs with fi — 10~^ can host no more than one or two equal-mass planetesimals per octave in radius, while cold discs with 
fi — 10^^ can host up to 12. 

6.2 Hot discs 

6.2.1 The maximum surface density for hot discs 

Comparing Figures[T]and|2]shows that a necessary condition for the survival of a hot disc is that the allowed regions in the two figures 
overlap, and that this requires (i) t'^ jv'^ < [i (the minimum planetesimal mass for which the collision time is longer than the age must be 
smaller than the disc mass) and (ii) t /u^^"^ < 1 (the minimum eccentricity for which the collision and gravitational scattering times are less 
than the age must be less than unity). The first of these can be written more accurately using equations l llOt and l l27t as 

M < Mmax.c (68) 

where 

7:'/'f'J' V 



/^max,c 



48/;V» (!ilo)>/= 



(69) 



To describe condition (ii) more accurately, we re-write the relaxation time i37b using equation 

= ( -T^\ ihQ-"^ + HQ"' + f.Q''' log A), (70) 



eo V^TrppM© 

in which we have assumed A ^ 1, as is usually the case, and neglected the contribution of shear-dominated encounters. The minimum of 
the expression in brackets occurs when O — /i ^'^/(S/s log A)^''^ for log A S> 1, and equals 0.97(log A)^/^. Since eo < /e (equation (8j), 
the relaxation time cannot be greater than the age to for any rms eccentricity unless 

/.i ^ /^max, relax (71) 

where 

„ , -44x10- rAV^V''7^V''7i^V''7^^V'7^VV^)' (72) 

Mmax.roiax - 4.4 X 10 g 5 J g gg J q 28 J ^10 7 \Z Z cm'^ ) \ZGyr) llO AU J ' ^ ^ 

The existence of a maximum surface density implies a maximum value for the IR excess emission due to dust (see i j8.2b . 



6.2.2 The maximum number of planetesimals in hot discs 



When a hot disc satisfies the constraints l |68b and ([TTJ, for a given surface density and semi-major axis there is a a minimum planetesimal 
mass and maximum number of planetesimals, given approximately by (cf. Fig.|2j 

mJSn _ (Snto)' ^hot _ /PpaN2 1 ^^^^ 



Mq MqpI ' V E y into 
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More precisely, by evaluating tc — to (O <^ 1) we get; 
hot _ 2304/1* (Snto)' 



Pp 



2WM, 



VlO-4/ 1^0.69/ V3gcm-V VSGyry VIOAU/ 



2304/f V E / {Q.to^'^ 



6.3 Warm discs 

Similarly, for given values of the surface density and semi-major axes there is a maximum mass im^™) - and thereby a minimum 
number {N^™)- of planetesimals in warm discs. Generally, the maximum mass and mininum number of planetesimals need to be evaluated 
numerically, but for the cases considered below ( i]6.4t , these extremes are attained when the line tc = to intersects the line S = 1 (cf. Figs. 
|4]and[5}. At this point. 




6.4 Sample discs 

As described at the end of [j2l we examine six possible planetesimal discs (TableO, with semi-major axes a = 1, 10 and 100 AU, and 
dimensionless masses ^ — 10"'' and lO^''. The allowed values of velocity dispersion Or and planetesimal radius r are shown in Figures|3] 
andg] In all cases, we assume that the disc age is to = 3 Gyr and the planetesimal density is pp = 3 g cm . 

A. a = 1 AU, = 10~*: Hot and warm discs cannot survive (see equations 1401 . 1681 , 1711 and \15V \. Cold discs can survive, but 
only for a narrow range of planetesimal masses (the thin blue trapezoid in the upper panel of Figure [5): there can be at most one or two 
planetesimals per octave of semi-major axis, of mass m™/^ ~ GQMq ~ 0.2Mjupiter (equation 1651 ). Such discs are rather similar to some 
of the many extrasolar planetary systems already detected by radial-velocity variations in the host star. 

B. a = 1 AU, /I — Hot and warm discs cannot survive. Cold discs can have planetesimal masses in the range O.l-lMg (equation 
1651 ): the lower limit corresponds to about a dozen planetesimals per octave of semi-major axis (equation 1671 ). These discs may be detectable 
with space-based transit surveys and are reminiscent of the terrestrial planets in our own solar system. 

C a — 10 AU, ^ = lO"**: Hot and warm discs cannot survive. As in the case of disc A, only one or two planetesimals (or planets) 
per octave of mass m'^l^ ~ 60 can survive in a cold disc. We suggest in i]7.5.3l that some discs of this type may be detectable by 
gravitational microlensing. 

D. a = 10 AU, jj, — 10^'': Hot, warm, and cold discs can all survive. The hot discs may contain up to about 5000 planetesimals per 
octave with minimum masses of m^f„ ^ 10^" g, about the mass of Ceres. The dynamical constraints on cold discs are the same as for disc 
B; such discs are not detectable with current or planned transit surveys because the probability of transits is too small and the orbital period 
is too large but could be detected by targeted searches for gravitational microlensing (see Figure|7j. A wide range of warm discs is possible, 
with at least 200 planetesimals per octave and masses at most 0.4 times that of the Moon. A more typical warm disc might have 5 x IC 
planetesimals per octave, of radius 100 km, with velocity dispersion ar — 0.1 km s~^ and a collision time of 0.6 Gyr. 

E. a = 100 AU, /i — 10""*: Hot, warm, and cold discs can all survive. As in the case of discs A and C, the cold discs can have only one 
or two planetesimals per octave. Hot discs can have up to 2 x 10^" objects per octave with masses of at least 4 x 10^® g. Warm discs contain 
at least 5 x 10* objects per octave with maximum masses of 10^^ g (r ~ 40 km; approaching the largest sizes of comet nuclei). 

F. a — 100 AU, n = 10~^': Hot discs have at most 2 x 10^^ objects per octave with m[^i'j ~ 4 x 10^'^ g. Warm discs contain objects 
with m < 10^^ g(r — 400 m), of which there are at least 5 x 10^^ per octave. The constraints on cold discs are the same as for discs B and 
D. 

It is worthwhile to emphasize that Figures [SJH] do not show how a planetesimal disc of a given initial mass will evolve, but rather 
whether a disc with a current mass Mdisc can survive in approximately its current state for 3 Gyr; the 'allowed' regions in these Figures can 
be interpreted as the allowed regions for discs with an age to — 3 Gyr on the assumption that it is unlikely to find objects in states that evolve 
on a timescale much less than their age. 
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Disc A: 0^1 AU ^^^10'" (M„3,^100 M,,,J, to^5 Gyr 




Figure 3. Allowed values of planetesimal velocity dispersion and radius for discs with age to = 3 Gyr, semi-major axis a = 1 AU, and dimensionless mass 
^ = 10~^ (A^disc — 100 M^) and fi = 10~® (A<fdisc — ^®), i-E-. discs A and B of Table[T] The blue shaded region denotes allowed cold discs. There 
are no allowed hot or warm discs. The various lines represent conditions for: gravitational stability (equations 1181 and 1231 ). gravitational scattering time 
exceeds disc age (equation 1311 ), thin disc (equation (s)), N > 1 (equation 1101 ). collisions are not erosive/disruptive (equations 1401 and 1441 ), and viscous 
spreading time exceeds disc age (equation 1431 ). The arrows attached to each line indicate the region in which long-lived discs could exist. The dividing line 
between hot/warm and cold discs is given by equation I14t and the dividing line between collision speeds greater than or less than the escape speed from the 
planetesimal surface (© < 1 or > 1 respectively) is given by equation {26). 
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Disc C: a=10 AU, /x^lO"" (M,i3,= 100 M,„,J, to=3 Gyr 



E 10* - 




r (cm) 

Disc D: Q=10 AU, ^JL=^Q-^ (M,,3,= 1 M,,,J, to=3 Gyr 




(cm) 



Figure 4. Same as Figure[3] but for a = 10 AU, i.e., discs C and D. For disc D, the pink sliaded regions denote allowed hot discs (solid color) and warm discs 
(vertical hatching). An additional constraint for disc D is tc > (equations 1271 and 1301 ). The dashed curve represents collision-limited discs (equation 

dii). 
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6.5 Collision-limited discs 

Collision-limited discs are shown by dashed curves in Figures [5] and |4] for discs D, E and F (equation 1571 ). The allowed ranges of 
masses are 10^* g < mmax < 8 x 10^* g (disc D), 4 x lO" g < m^ax < 3 x 10^^ g (disc E) and 4 x 10^^ g < mmax < 4 x 10^** g (disc 
F). 



7 DETECTION TECHNIQUES 
7.1 Radial velocity measurements 

The stellar wobble or reflex radial velocity induced by an edge-on disc containing planetesimals of mass m is approximately 



''^wobble 



GMqN m 
a Mn 
fn 



GMq (76) 
aN 



= 3.0 m s"^ /„iV" 



1/2 f a A'^/y M 



10 AU/ VIO 



-1-4 



With current technology we can detect reflex velocities as small as «wobbic ~ 1 m s ^ with orbital periods as long as ~ 10 yr, 
corresponding to a ~ 4.6 AU. For hot discs with ages of several Gyr, equations l l71t and ( 176b imply t^wobbic < 0.1ms-^(a/10AU)^''^too 
small to be detectable. 

Cold discs are detectable if the semi-major axis is small: combining equation i67b with the second of the equations above, we have 

^wobble > 2.3 m s U (yoAuJ llF^j • ^^^^ 
Thus all long-lived cold discs with fi > 10~* and a < 10 AU can be detected by current radial-velocity surveys; this of course is because 
gravitational stability requires that they have only a few large planets. 



7.2 Transits 

Space-based transit surveys such as NASA's Kepler mission are capable of detecting photometric variations as small as ~ 10~^, 
corresponding to the transit of a planetesimal of radius 2000 km — 0.35-R®. Thus edge-on discs containing lunar-mass planetesimals may 
be detectable by transit surveys. Reliable detection of transits of such small objects requires that the stellar variability is negligible; this is 
likely to be true for at least some stars as the solar variability on the hourly timescales relevant to transit detection is only a few times 10^^ 
i Batalha et al]|2002h . 



Other criteria for detectability of edge-on planetesimal discs by transits include the following: (i) The orbital period must be less than a 
year or so, so that several transits of a given object can be detected in a mission of reasonable duration, (ii) There must not be too few transits, 
that is, at least one planetesimal in the disc must transit the star. If the disc is nearly edge-on and the characteristic thickness h (equation 
t4\ ) is small compared to the stellar radius i?* then most planetesimals transit the stellar disc in the course of an orbit, while if ft ^ _R* the 
fraction of transiting planetesimals is (2/11)^ ^'^ R* /h, so the expected number of transiting planetesimals is roughly 



Nt = N min 



(78) 



(iii) There must not be too many transits: if multiple planetesimals are transiting the disc at any one time, the fluctuations in stellar flux will 
be difficult to distinguish from normal stellar variability. The average number of planetesimals in transit at a given time is NRi,/{na) if 
h <ti R* and NR^ / (y/Sivha) if Ri, so the average number of planetesimals in transit at any instant is 

nt=Nmm( — ,^^]. (79) 



7ra WgTTha 



7.3 Microlensing 

7.3.1 Microlensing by individual planetesimals 

The classical lensing equation is (e.g., ISchneider et alj|l992l) 

adi, — 



d L ' 

where Rs — 2Gm/(? is the Schwarzschild radius of the lens and L is the projected separation between the light ray and the lens in the lens 
plane. The distance to the source and lens and the separation between them are denoted by d*, d and di* = d* — d, respectively. The quantity 
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a is tlie subtended angle between the lines of sight to the lens and the source. The Einstein radius is defined by the value of L when a = 
(i.e., lens and source are aligned). 



-\/Gmd*C(l-C) 



(81) 



where C, = d/d^. 

A planetesimal can be far from or near to its parent star, where 'far' and 'near' are defined with respect to the stellar Einstein radius. 



Re = 4.0 AU 

with Mt being the stellar mass. 

The magnification of the total flux from the source is 



rf. C(1~C) 



-,1/2 



MoSkpc 0.25 



A = 



+ 2 



(82) 



(83) 



itVu^ + 4 ' 

where u = ad/vE- In our simple treatment, we assume that a planetesimal can produce a detectable lensing event of non-negligible magni- 
fication when the distance ad between the lines of sight to the planetesimal and the source star, measured in the lens plane, is less than the 
Einstein radius; this corresponds to u = 1 or amplification A — 3/^5 ~ 1.34. 

We must check that the planetesimal radius is small compared to the Einstein radius to ensure that the magnified light curve is not 
blocked ( Agol 2002). We have 

1/2 



r 
rE 



3c" 



16nGpprd-,C{l - () 



= 0.66 



d. C(l-C) 



lkm3gcm-3 8kpc 0.25 



-1/2 



(84) 



so obscuration by the planetesimals is unimportant if they are much larger than 1 km. 

In addition, we require that the stellar (source) radius as projected on the lens plane - equal to the stellar radius i?* multiplied by ■ 
cannot be much larger than the Einstein radius, to ensure that the magnified light curve is not smeared out. We have 
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(85) 



We have parametrized in terms of the solar radius since this is the typical size of the source star in existing planetary microlensing events 
(at 8kpc the corresponding angular size is 0.6 ^as). With the nomina l parameters and ( — 0.5, ii*C/^E < 1 only for planetesimal radii 
r > 3750 km, corresponding to mass m > O.llM® jPaczvnski|[l99i gives a similar estimate, O.O7A/0). This limit is conservative because 
the magnification of extended sources remains substantial when the projected source radius is as large as several times the Einstein radius 
- for example, a unifor m source whose centre is separa ted from the lens by one Einstein radius is magnified more than a point source so 
long as R^^C/rE < 2.17 dCouldf 1994l : lwtt & Maolll994h . Thus, microlensing sea rches are likely to be sensitive to planetesimals as small as 
~ W'^-''M(s or a few times the mass of the Moon (but see' Heng & Keetonll2009lFI . 

If the transverse velocity of the lens relative to the source is v± , the characteristic duration of the event is 

1/2 



2rE 

v± 



16 min 



Vj_ 



100 km s 



3/2 



Pp 



3gcm-='8kpc 0.25 



(86) 



.1000 km/ 

In most cases, the transverse velocity is dominated by the apparent angular speed of the source star relative to the host star of the planetesimal, 
rather than the motion of the planetesimal around its host star. 

If we assume that the surface density of the disc is uniform over a circle of radius a, the probability of lensing at any given moment for 
a star whose image lies within the disc (i.e., the optical depth) is 

A* d* 



10 AU 



C(i-C) 

0.25 ' 



(87) 



^><10"^^"'VlO-8kpc 
The optical depth of a planetesimal disc of a given size a and mass Mdr. 

mals. Thus (for example) the optical depth for a disc composed of 100 Earth-mass planets (Afdisc/Af© = 3 x 10" 
as the optical depth of a single 0.3 Jupiter-mass planet at the same radius. 



7r/mSa is independent of the mass of the individual planetesi- 

, p — 10~*) is the same 



^ Events associated with lower amplifications have lai'ger microlensing cross sections, i.e., ■k{4>'''e)'^ where </> > 1. Planetesimals with Einstein radii smaller 
than the projected size of the source star may contribute appreciably to the expected number of events per planetesimal disc crossing, because the range of 
masses involved in microlensing now extends down to much lower values. 
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A related quantity is the probability that at least one lensing event by a planetesimal will be seen at some time during the passage of the 
source star near the host star of the planetesimal disc. If the impact parameter associated with this passage is small compared to the size of 
the planetesimal disc, this probability is 1 — exp(— Mens), where 



Nlens = 



8 



TV 



1 1/2 



C) 



(88) 



= 9.1 X 10" 
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Notice that the probability of observing an event goes up as the planetesimal mass goes down, since A'^ ~ Ea^/m (although the duration of 
the event is shorter). In this respect, a disc consisting of many small planetesimals may actually be easier to detect than a single large planet. 

The short duration | |86I ) of the events is one of the principal challenges in reliably observing microlensing by planetesimals of an Earth 
mass or less. To avoid being swamped by noise it is useful to focus on source stars that are experiencing - or have recently experienced - 
microlensing by an intervening star. Strong amplification by the host star of a planetesimal disc requires that the impact parameter is less than 
the Einstein radius of the host star, given by equation ( |82b . Since many planetesimal discs may be substantially larger than Re the source star 
should be monitored for short-duration events for some time after the amplification by the host star has returned to unity. 



7.3.2 Other microlensing effects 

Planetesimal discs can produce other signals in microlensing surveys. Zheng & Menard i 2005h point out that if the source star for a 
microlensing event hosts a debris disc, the mid/far-IR light curve will contain a component determined by the surface-brightness profile of 
the thermal emission from the debris disc. Similarly, the optical and near-IR light curve will contain a component from the scattered light 
from the disc. In both cases we may expect that the light curve is no longer wavelength-independent. These effects are challenging to detect 
since (i) many debris discs are much larger than the stellar Einstein radius Re — 4 AU (equation 1821 ). so the maximum magnification is 
only ~ {Re/o,)'^', (ii) accurate mid/far-IR photometry is exceedingly difficult, except from space; (iii) the fractional flux of scattered light is 
small, typically 10~^ to 10~^ in observed debris discs (e.g.,[wvatt 2008). 

Other signals may arise if the lens star hosts a planetesimal disc. The overall mass distribution in the disc will contribute to the magnifica- 
tion and thereby distort the microlensing light curve, but this distortion will be difficult to detect because the magnification due to the disc will 
only be of order /x7?b /a where /i is defined in equation Q; we have assumed tha t the disc is not far from face-on and that its semi-major axis 
a is larger than the Einstein radius of the host star (see Hundertmark et all200 9 for a discussion of lensing by edge-on discs). A potentially 
more sensitive probe is high-magnification events. The gravitational field from distant stellar companions or ot her external mass distribu- 
tions can produce a characteristic double-peak structure near the point of maximum magnification; for example, Kim et al. 1 200 sh estimate 



that typical high-magnification events can detect stellar companions with mass m and separation Ad such that tti/Mq > (Ad/100 AU) . 
Unfortunately, these events are insensitive to dist ant discs in most case s, for the following reason: the deflection angle from a surface mass 
density distribution Eproj (x) on the sky plane is i Schneider et al.|[l992 ) 



4GSp,oj(x') x-x' 



(89) 



which is proportional to the gravitational field from a cylindrical mass distribution with density p(x\,X2,Xi) oc Eproj (a;i, 2:2). The deflection 
angle for light rays passing inside an inclined ring of material is therefore proportional to the gravitational field inside an elliptical, cylindrical 
shell, which is zero from Newton's theorem. Only discs in which Eproj is non-zero near the host star (e.g., nearly edge-on discs with a 
significant thickness) will affect the light curve near the peak magnification. 



7.4 Infrared emission due to dust generated from collisions 



Most extrasolar detections of planetesimal discs are based on measurements of IR excesses, i.e., the presence of a debris disc jWvatt 

I2OO8') in which a steady supply of dust is generated by planetesimal collisions; the dust is heated by the host star and the heated dust generates 
IR emission. At a given wavelength, the IR excess, /ir, is the dust luminosity measured relative to the stellar luminositjljl In this subsection 
we estimate the IR excess due to planetesimal collisions in discs. The discs considered here differ from the collision-limited discs of ^in 
that the collision time exceeds the age of the discs, tc > to- 

If the dust grains are treated as gray bodies, their equilibrium temperature is 

IbTra^crsB/ V " / V-^© 



Tdust 



(90) 



^ Confusingly, / is also used by many authors to denote the bolometric luminosity of the disc relative to the stai'. For particles with absorption efficiency 
Qa = 1 (equation 1921 ). this ratio is equal to the geometrical optical depth Tp defined in equation I12t . 
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Figure 6. Allowed values of the infrared excess, /ir(A), for hot discs around solar-type stars with masses M^^^c = ™d 10 A/®. The shaded regions 
denote hot discs that are detectable via their IR excesses; the assumed detection thresholds /dot at A = 24 and 70 im\ are shown as horizontal lines. 



where ctsb is the Stefan-Boltzmann constant. The peak wavelength of the black-body spectrum \B\(\, T) at this temperature is Amax = 
13 Mm (a/ AU)^/^(£/£0)"^/*. The IR excess is then 

where and i?* are the stellar temperature and radius, Qa is the absorption efficiency, and dN^ust (r) is the number of dust particles as a 
function of their radius r. As discussed at the end of Sj5] the absorption efficiency can be approximated as 

=min{l,X} where X = (92) 

A 

so if the number of dust particles is a power law in radius, 

diVdu.t = Kr-'dr, (93) 

we have 

The radius exponent q is related to the mass exponent p defined in equation ( 145 l l by g = 3p — 2. For small particles such as dust, we 
expect that the specific kinetic energy required for disruption, Qjj, is independent of mass (cf. i]4.1t . In this case p — 11/6 (equation l lSOt i 
so g = 7/2, and we shall use this value in evaluating equation l |94| ( numerically. Note that equation ([94} is only valid if 3 ^ g ^ 4 and if the 
minimum grain size in the distribution is much smaller than A/27r. 

The physical processes governing the effects of collisions are outlined in i]4.II In a steady state, the rate of dust mass production in a 
monodisperse planetesimal disc of the kind we are considering is 

$™ - ^, (95) 

where as usual m and N are the mass and total number of planetesimals and tc is the collision time, given by equation ( I27t . To estimate the 
corresponding dust mass, we use the conservation of mass flux. 

First we generalize the collision time Ml\ to the case where particles of radius ri are colliding with particles of radius r2. We have 

(ri) = 2 ' TT ' /iCTr / (ri+r2) — - dr2; (96) 

here we have assumed that self-gravity is negligible {Q <C 1). If the number density is a power law in radius, dno{r)/dr oc r"', and the 
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integral is dominated by projectiles with radii r2 much less than the target radius n, we have 

t-(n) = ^^^^a..?.-^ (97) 
q — 1 dri 

where 6 = ri/rmin ^ 1, i.e., ri S> rmin- If we define rmin to be the minimum projectile radius that will disrupt a grain of radius n, then 
the mass flux in the disc is roughly 



where m(ri) — |7rpprf is the mass of a grain of radius ri. The number density and the total number of grains are related by equations ^ 
and 

fdN^\ ^j^^.~,^^./2^mf^}^^^^^:Mr^^ (99) 
\d\ogr J eo il dri 

Equating equations l |95l l and l |98t . using equation l |99t to eliminate dno{ri) /dri in favor of K, and assuming that the velocity dispersion, 
internal density, and radial distribution of the dust and planetesimals are equal, we have 

e'^'^^K'^ « ioivV^(i + f2<d/h), (100) 

which relates the properties of the dust distribution on the left side of the equation to those of the parent planetesimals on the right. Although 
we have assumed q — 7/2 for the dust particles (i.e., the particles that dominate the IR emission), which is equivalent to the assumption 
that the specific dispersion energy QJ, is independent of mass for these particles, this derivation does not require any assumption about 
the dependence of Q*^ on mass for larger bodies - the argument relates the mass flux from the destruction of equal-mass planetesimals 
(equation 1951 ) to the mass flux in dust (equation 1981 ) and since mass flux is conserved the properties of intermediate-mass bodies are 
irrelevant. The value of 9 is estimated from the specific dispersion energy QJ,; following the discussion after equation i44\ a typical value is 
Qo ~ 10^ erg g~^ and we equate m(ri)(3o to im{rmin)o-r to obtain 



„2/3 



-1/3 



(101) 



10 AU 10^ ergg-i^ 

The approximations that we have made to derive equation l IlOOI l are invalid unless 9^1. 

Equation dlOOl ) determines the normalization K of the dust distribution, which is substituted into equation l |94b to determine the IR 
excess. This derivation is for hot discs, and implicitly assumes a monodisperse disc in which collisions slowly feed a population of smaller 
debris. These assumptions are only valid for discs in which the collision time exceeds the age, tc > to. Once tc ~ to the appropriate model 
is a collision-limited disc ( ij4.3t , and the factor tc in equation l |95t should be replaced by to. The analogue to equation l IlOOI l is then 

32Q,to fi pp 

As discussed in i]3.2.2l collisions also occur in warm discs, but in this case the collisions do not disrupt the planetesimals and a collisional 
cascade is not established. We evaluate the IR emission properties of warm discs in i]7.5l below by assuming that most of the emission comes 
from the planetesimals themselves. 

For calibration-limited detections - in which the limiting factor is the accuracy of the extrapolation of the photospheric flux to long 
wavelength^ ~ the minimum detectable flux produced by th e dust, normalize d by the stellar flux, is /dct- For illustration, we set f^ct = 0.1 
(A = 24 fim) and 0.55 (A = 70 /im), similar to the limits in ISu et all lioodi) . We show examples of hot planetesimal discs with detectable 
IR excesses in Figure[6] It is apparent that the IR excess is a poor diagnostic for the disc mass A/disc : the infrared flux from discs of a given 
mass and semi-major axis in Figure|6]can vary by more than an order of magnitude. 



7.5 Probing disc mass and size 

We now ask what long-lived planetesimal discs are detectable by the methods we have discussed in i) i]7.1tf01 To efficiently explore the 
four-parameter space of disc mass and semi-major axis, planetesimal radius, and velocity dispersion (Afdisc, a, r, Ur), we randomly generate 
3 X 10" discs, uniformly sampled on logarithmic scales: 10~* ^ Mdisc/M® lO", 0.1 a/ AU sC 1000, 10"" sC r/cm s; 10^^ and 
cm s ^ ^ 10^. We then ask whether each disc can survive for to = 3 Gyr and is detectable by one or more methods using the 
detection thresholds described below. Figure [7] shows the detectable planetesimal discs as projected onto the Mdiac-i plan^. 



« See ^2.4 of lwvatj j2008l) for a discussion of calibration- vs. sensitivity-limited surveys. 

^ Note that the density of generated points in Figure|7]is generally lower at higher disc masses (Mdisc It 100 Af ^ ), which is a surprising result since massive 
discs should be easier to detect. The low density comes about because the range of allowed planetesimal radii for cold discs becomes narrower for higher disc 
masses (see i|6.4l and Figures[3]|4]and[5j, and we are sampling logj^o '" unifoiTnly. 
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Figure 7. Different detection tecliniques probe different ranges of planetesimal disc mass and semi-major axis. Eacli of the four disc parameters (mass, semi- 
major axis, planetesimal radius, and velocity dispersion) is randomly generated and the various detection criteria are checked (see text). Only discs that survive 
for 3 Gyr are shown. The detectability criteria are summarized in ii il7.5.Ht7.5.5l 

7.5.1 Radial velocities 

We consider a planetesimal disc to be detectable by this method if the orbital period 2n{a'^ /GMq)^^^ is less than ten years and 
■Uwobbio 5S Vdct, where Uwobbic is given by equation ( 176b and the detection threshold Vdct = 1 m s~^. Black crosses in Figure|7]denote discs 
that survive for 3 Gyr and are detectable by this method. The minimum detectable mass scales oc a^^^, as expected. No warm or hot discs 
of age 3 Gyr were detectable by this method. Not surprisingly, planetesimal discs that are detectable by radial velocity variations in the host 
star tend to be massive and contain a small number of large bodies, i.e., planets; they resemble Disc A of %.4\ 

7.5.2 Transits 

As described in i]7.2l planetesimals in nearly edge-on discs can be detected transiting their parent star if r > 0.35R^, and a < 1 AU. 
We also require that Nt > 1 (equation 1781 ) and nt < 1 (equation 1791 ). It could be argued that the condition A'^t > 1 is unnecessarily 
stringent, since even if Nt 1 a. fraction of discs with these properties could be detected in a large transit survey. 

No warm or hot discs of age 3 Gyr were detectable via transits. Discs A and B of i]6.4l are detectable via transits. 

7.5.3 Microlensing 

We consider a planetesimal disc to be detectable via microlensing if the microlensing optical depth ricns (equation 1871 ) exceeds 10~^ 
(for comparison, the measured microlensing optical depth towards the Galactic bulge is 2-3 x 10"®). Since ricns is maximized when the 
planetesimal is halfway to the source, we adopt C, = d/di, = 0.5 for illustration. We also assume a solar-type source star at a distance 
= 8 kpc. With these parameters, no warm or hot discs are detectable by microlensing. In Figure[7l the cutoff for Afdisc ^ O.IM^ arises 
because the projected source size becomes larger than the Einstein radius, while for a > 1 AU we have the detectable disc mass scaling cx a? 
(equation 1871 ). Discs A, B, and C of i]6.4l are detectable by microlensing. 

7.5.4 Debris discs 

Debris discs are dynamically hot discs that produce a collisional cascade of dust whose associated IR excess exceeds the detection 
threshold, /det- We consider both hot and collision-limited discs. As discussed in i)7.4l we take /dct =0.1 (24 ^m) and 0.55 (70 /im). 
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In Figure [T] there is a minimum and maximum detectable disc mass at a given semi-major axis. Tiiis feature appears to arise because tlie 
relaxation time trdax (equation IB?! ) at fixed disc mass and velocity dispersion has a minimum near 8 = 1. There is also a cutoff at small 
semi-major axes, which arises because the maximum allowable surface density for hot discs (equations 1681 and 1711 ) is a strongly increasing 
function of semi-major axis. Discs E and F of ii6.4l are detectable via their IR excesses at 70 /im. 

It is remarkable that the discs detectable by radial-velocity/transit surv eys or micr olensing do not overlap with those detectable from 
IR excess. This result is consi stent with the obser vational findings of Beic hman et al. (2005), Bryden e t al. (2006), Greaves et al. (20(3), 
2007 ), and Kospal et alj ( I2OO9I) that there is little or no correlation between the occurrence of planets and debris discs. 



Moro-Martm et al 



This lack of overlap does not preclude the possibility that a single host star may have planetesimals that are detectab le by both methods, 
so long as the planetesimal disc extends over several octaves in semi-major axis. Beichman et al. (20o3) and Kospal et al. ( 2009) list six 
and ten pla net-bearing stars w ith debris discs, and planets have been imaged in the debris-disc systems HR 8799 jMarois et aljboosi) and 
Fomalhaut fealas et alj|2008l) . 



7.5.5 Warm discs 

In warm discs collisions may cause cratering of the planetesimals but do not shatter them. Although cratering collisions produce signif- 
icant amounts of dust, they do not establish a coUisional cascade of the kind described in i)4.2l so it is likely that IR emission is dominated 
by the planetesimals themselves. Based on this assumption we show the detectability of warm discs from their IR emission in Figure|7] We 
again adopt /dot = 0.55 at 70 /im (we choose not to show detectable warm discs at 24 /im so as not to over-crowd Figure|7](. 

In the Figure we see that there is again a minimum and maximum detectable disc mass at a given semi-major axis. This feature is now 
associated with the constraints set by non-gravitational forces (tpR > to and /3 < 0.5; see [j5} and condition ( 1401 ) for warm discs, respectively. 

In principle, cold discs may also be detectable through the IR emission from the planetesimals, but we found no such discs given our 
assumed detection limits at 24 /im and 70 /im. 

7.5.6 Undetectable discs 

Many planetesimal discs that survive for 3 Gyr are not detectable using any of the methods described in this section. 



8 DISCUSSION AND SUMMARY 



8.1 Can warm discs mimic debris discs? 

It is generally believed that the IR excesses around main-sequence stars older than a few Myr are due to dust that is produced in 
collisions between large solid bodies orbiting the star (hence the term 'debris disc'). Direct evidence that the emitting material is dust comes 
from several sources: (i) Submillimeter observations of a handful of debris discs show that the absorption efficiency Q g (equation 1911 ) 
declines roughly as for wavelengths > 100 tim ' Dent et alj|200ol : Williams & Andrews 20061 : Backman et aljEoOT ). suggesting grain 
sizes of a few tens of /im (equation 1921 ). (ii) Chen et aLn200^ ) obtained Spitzer Space Telescope infrared spectra of 59 stars with IR excesses 
and found five with 10-20 /jm features that imply the presence of micron-sized silicate grains, (iii) The pol arization of scattered light from the 
debri s discs around /3 Pictoris and AU Microscopii is consistent with simple models of scattering by dust ( Gledhill et aljl991 : Graham et al 
2007h . 



Despite this evidence, it is instructive to consider the possibility that in some stars the IR excess arises not from dust produced by a 
coUisional cascade but rather from a population of planetesimals with much larger radii. The most likely candidates are warm planetesimal 
discs, in which the collision time is less than the disc age but the collision velocities are too small to destroy the planetesimals over the 
lifetime of the disc. To simplify the calculations, we consider the lowest possible radial velocity dispersion for warm discs, which occurs 
when the Safronov number = 1 (cf. discs D, E, and F in Figures |4]and[5l(. Using equations ( I25t . ( 1271 ) and the third condition in equation 

we obtain 

1 ( Qh 



r s; 



(103) 



8(/i + /2) \GT.Q.to, 

If 6 < 1, then the preceding constraint becomes stronger, i.e., the numerical coefficient in equation ( I103t becomes larger. 

The minimum planetesimal size in warm discs is typically set by Poynting-Robertson drag. Using equations ( 158b and ( I63t , the condition 
ipR to yields 



r > 0.7cm p0 



to 



3 Gyr 



Pp 



(104) 



^3gcm-V VlOAU/ 

In Figure[8]we show the constraints JI03t and JI04t for a = 1, 10, and 100 AU; also shown are lines of constant optical depth Tp 
(equation II2I ). This optical depth is equal to the ratio of the bolometric disc luminosity to the bolome tric stellar lu minosity and hence 
provides a convenient measure of the detectability of the disc. Known discs typically have Tp > 10~^ ( IWvatlll2008l) . We conclude from 
Figure [8] that the IR emission from so-called 'debris discs' at a = 100 AU could in some cases be coming from planetesimals as large as 
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Figure 8. Allowed surface density of warm discs as a function of the planetesimal radius according to equations )103t and )104t . The allowed region is shaded 
in different colours for different values of the semi-major axis a = 1, 10 and 100 AU. The disc age is assumed to be 3 Gyr. Also shown are lines of constant 
optical depth Tp (equation 1121 ). 



r ~ 10 m (m ~ 10^ g). A strong test of this possibility is that the emission spectrum from such a disc should resemble a black-body 
spectrum, even at submm wavelengths (or a superposition of black-body spectra if the emission originates from a range of disc semi-major 
axes). 



8.2 The maximum optical deptti of a debris disc 



Wvatt et al.l ( 120070 argue that a simple model for the collisional evolution of planetesimal discs implies that the maximum optical depth 



or fractional bolometric disc luminosity is (their equation [21]) 

) = 1.6 X 10-^ ( — ) (105) 

This result is based on several plausible but arbitrary assumptions (planetesimal radius r = 2000 km; rms eccentricity eo = 0.05, strength 
Qi, = 2 X 10^ ergg'\etc).To examine the applicability of this result, we employ the same Monte Carlo procedure used to produce Figure 
|7]to randomly generate hot, collision-limited, and warm planetesimal discs and calculate the corresponding values of Tp. For warm discs, the 
optical depth is given by equation J12l ( with A'^ and r equal to the number and radius of the planetesimals. For hot and collision-limited discs, 
the optical depth is given by 

Tp = Qar'dAdust; (106) 

taking the absorption efficiency Qa from equation J92t and the number of particles dN^ust from equation l |93l l with q — 7/2, we find 

K f -'/^ 



^^ = ^UJ ' ''''' 
where K is taken from equations JlOOt or ( 1102 1 ) for hot and collision-limited discs, respectively. In the discussion below we assume A — 
70 urn. 

The results are shown in Figure|9]for disc ages to — 30 Myr, 300 Myr and 3 Gyr. We also show Wyatt et al.'s estimate dlOSt for to = 30 
Myr as a dashed white line; this is easily scaled to other ages since rp™^"' oc l/to- Wyatt et al. estimate the uncertainties involved to span 
~ 2 orders of magnitude and this is reflected in the light blue band shown in Figure[9] 

We are able to generate hot and collision-limited discs with optical depths substantially larger than the estimate of equation U05t ; 
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Figure 9. Optical depth associated with dust grains in hot (left) and collision-limited (middle) planetesimal discs, as well as for planetesimals in warm discs 
(right). Blue, green and yellow symbols are for planetesimal disc ages of 30 Myr, 300 Myr and 3 Gyr The white dashed line and light blue band represent 
the maximum optical depth and its associated uncertainty, respectively, as estimated by Wyatt et al. (2007) - see equation U05t . Also shown are the obsei-ved 
values for HD 12039 and BD +20 307. 



5 rou gl 



axis is quite different. The numerical simulations of lLohne et alj ( 120081) also show that the scaling of Tp with a is generally more complicated 



however, these still lie within the estimated range of uncertainty given by IWvatt et alj j2007l) . A major source of uncertainty is in the 
planetesimal strength Q*q. In the calculations shown in Figure[9]we used a mass-dependent Q n as defined in equation l l44b . As a check, we 
carried out additional calculations assuming the constant value adopted by Wvatt et ^ 1 2007 ) (QJ, = 2 x l(f erg g~^), and found that in 
this case our results (not shown) agreed more closely with equation ilOSb . 

Our simulations show that the maximum optic al depth is roughly e x t^^,as, predicted by equation ( 1105b . but the scaling with semi-major 
is quite different. The numerical simulations of lLohn 
than a power law (see top right panel of their Figure 11). 

We also find that the maximum optical depth of warm discs can be almost an order of magnitude higher than that of hot and collision- 
limited discs of the same age and semi-major axis. The absence of warm discs in the bottom right corner of Figure|9]is simply a consequence 
of the defining condition of warm discs, tc < to (equation 1271 ). together with equation ( 164b relating the collision time to the optical depth. 
Discs exist below this line, but we label them 'hot' rather than 'warm'. For collision-limited discs, there is a similar cut-off caused by the 
thin- disc condition (equ ation Isl ) and equation ( 157b imposing a maximum value for mmax (e.g., see Figure[5](. 



Wvatt et al.l ( 120070 point out that a nu mber of debris discs with a < IQ AU have optical depths that exceed the limit ( 1105b by factors 
of 10^ or more (see also Moor et al. 2009h . In Figure |9l we show tw o debris discs with small semi-major axes, HD 12039 {to = 30 Myr) 
and BD +20 307 (to = 300 Myr), taken from Table 1 of IWvatt et al.! ( |2007|) . These stars have a range of spectral types, from F2 to K4, but 
our models based on a solar-type host star should still be reasonably accurate. We verify that the optical depth of BD -1-20 307 exceeds the 
maximum allowed for steady-state hot and collision-limited planetesimal discs with an age equal to the stellar age, while the optical depth of 
the disc around HD 12039 is consistent with steady-state models. 

In Figure [Tol we show all seven debris discs listed in Table 1 of lWvatt et aLl ( l2007l) . For each system, we use the quoted values of the age 
to and semi-major axis a to compute the maximum value of Tp for hot, collision-limited and warm discs. The Figure shows that two systems 
(HD 1 13766 and HD 12039) have optical depths consistent with a steady-state hot or collision-limited disc; one (BD -1-20 307) has an optical 
depth that is inconsistent with a steady-state hot, collision-limited, or warm disc (by factors of 100, 50 and 10, respectively); and four (HD 
72095, HD 69830, i] Corvi, and HD 98800) are consistent with warm discs but not hot or collision-limited discs. However, of these last four, 
the first three have 10 fim silicate features in their spectra which imply that the IR emission comes from micron-sized grains, thus ruling out 
warm discs as well. All of our conclusions about hot and collision-limited discs are consistent with Wvatt et al] i 2007 ), who suggest that the 
dust arises from planetesimals that have been scattered to small semi-major axes from a disc at much larger radii. 
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Figure 10. Maximum optical deptlis for hot and warm discs as calculated by our model, compared to the observed values, for the 7 systems listed in Table 1 
ofWyattetal. (2007). 



8.3 Summary 

We have described a unified model of the evolution of gas-poor planetesimal discs, which is general enough to apply to all Keplerian 
discs of solid bodies, including debris discs, asteroid belts, and planetary systems. Our model includes such processes as gravitational stability, 
evolution due to dynamical chaos, gravitational scattering, radiation and stellar wind pressure, Poynting-Robertson drag, and erosion or 
destruction by physical collisions. We characterize the discs by four parameters: disc mass (Mdisc), disc semi-major axis (a), planetesimal 
size (r) and radial velocity dispersion or rms eccentricity (ar or eo). The salient conclusions of our study include the following: 

• Planetesimal discs can be categorized as dynamically 'hot', 'warm' or 'cold' depending on whether the planetesimal orbits cross and 
therefore collide and whether the collisions are erosive/disruptive. In cold discs the orbits do not cross and collisions do not occur; in hot 
discs the orbits cross but the collision time is longer than the disc age, and in warm discs the collisions are frequent but gentle enough that 
they do not substantially erode the particles within the age of the disc. 

• Massive discs with small semi-major axes can only survive for Gyr timescales if they are cold. For example, after 3 Gyr hot discs at 
1 AU or 10 AU cannot exceed 1.3 x 10~*A/® or l.SM® respectively (see i]6.2b . Gravitational stability imposes an upper limit on the number 
of planetesimals per octave that can be present in a cold disc of given surface density; for example, a cold disc of mass lOOM® cannot host 
more than 1-2 planetesimals per octave, while a disc of mass IM® can host ~ 10 per octave (equation 1671 ). 

• Warm discs can survive for Gyr timescales over a wide range of semi-major axes and masses. At 1 AU warm discs that survive for 3 
Gyr must have mass < 10"'' M®; in this case the planetesimal radius is only 1 m, and warm discs composed of larger planetesimals must 
have even smaller masses (Figure [8}. At larger semi-major axes the allowed masses of warm discs and the planetesimals within them are 
much larger (Figure some cases warm discs may be detectable from the IR emission from the planetesimals themselves. 

• Planetesimal discs can be detected by a wide variety of observational techniques, including transits, gravitational microlensing, radial- 
velocity variations, and 'excess' IR emission ('debris discs'). With current technology the discs that can be detected by any of the first three 
methods are disjoint from those that can be detected in the IR (see Figure[7](. Many possible long-lived planetesimal discs cannot be detected 
by any method at present. 

Despite the length of this paper, our analysis suffers from several shortcomings. The assumption of a monodisperse planetesimal disc is 
oversimplified, and probably incorrect given our limited understanding of disc formation. We suspect that our results are reasonably accurate 
provided that the total mass in the disc is dominated by planetesimals in a relatively small mass range, but this suspicion should be tested 
by analysis of discs with a range of planetesimal sizes. Our results also depend on a number of poorly determined parameters of order 
unity (Table |2ll and do not incorporate a realistic model of the radial structure of the disc. In this paper we have deliberately ignored all 
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considerations of the formation process of planetesimal discs. It remains to be determined, by observations and theory, which of the wide 
variety of possible long-lived planetesimal discs are actually found in nature. 
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APPENDIX A: MODIFYING VARIOUS FORMULAE IN DONES & TREMAINE (1993) 



The rate of mass accretion in a rotating disc of planetesimals has been evaluated by iGreenzweig & Lissauen ( 119921) and Dones & 
Tremaine (1993; hereafter DT). We need to modify their formulae, because their results are for the mass accretion rate of a large body 
(a planet) on a circular orbit in the midplane of the planetesimal disc, while we are interested in the rate for a typical planetesimal in a 
monodisperse planetesimal disc. Unless otherwise mentioned, the notation used in this Appendix is the same as in the main text. 

In the dispersion-dominated regime, the rate of mass accretion is given by equations (72) and (90) of DT: 

r2.7603Eni?^, e«l, ^^^^ 

\6.0828E!:^^i^pi^^a-^ e>i, 

where Rp is the radius of the planet, i?H = a{M /MqY^^ is its Hill radius (as defined by DT, which is different from the definition in the 
present paper) and Gt is the planetesimal velocity dispersion in the radial direction. In the shear-dominated regime, we use equations (83) 
and (75) of DT: 



' ly^RT^ <Tr<n^/R^. 



We first write these formulae in terms of the number density of planetesimals in the midplane, no, using E = \/2TTnomaz/i^ (see text 
below equation [18] of DT): 



M : 



2.7603V27rnomcr^7?;, 



6. 0828V 2Tvnoma RpRl^a; , 

I — 2 —1 (^^3) 

10.lv 27rnom(Jzr2i?p_Rjj'^r , 

, 6.47V2^nom(J^i?p''^i^|''^ 
The collision time as defined in the present paper is = M /m. The following modifications are made: 

Rp 2r, (A4) 

RH = a {M/Mq)^'^ a {2m/MQ)^^^ . 

The first modification comes from assuming the colliding bodies have the same velocity dispersion, as opposed to one of them being on a 
circular orbit. The second and third modifications arise both colliding bodies have the same radius rand mass m, as opposed to one large 
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(A5) 



body having radius R and mass M wiiile tlie other has negligible mass and radius. We also note that CTz = a,, (io/eo ) , where we again choose 
io/eo = 0.5. Thus, the reciprocal of the collision time is 

'2.7603 X 4 X ^nor'^ar, 
6.0828 X 2 X ^noVi^ra^a-^ (m/MQ) , 
10.1 X 2^/® X ^noQra'^ (m/Me)^''^ , 
6A7 X 2 X 0Fnor^''2„3/2^^ (m/Me)^/^ . 

Finally, we decrease by \/2 since the vertical motions of the particles imply that the mean density is reduced by this factor compared 
to the midplane density. We also replace no by ^2/TvAffi/a-r: 

'2.7603 X 4 X Afflr^, 
6.0828 X 4 X AfQr^e, 
10.1 X 2^/** X AAfiVaV-i {m/M@f^^ , 



tc = < 



(A6) 



6.47 X 2 X A/'fir^/^a^''^ (m/M©)^''^ . 



Comparing with equation l|27ll, we get /i = 4 x 2.7603/16 = 0.690 and /2 
equation Oil, we get = 2^/'' x 10.1 = 22.67 and /s = 2 x 6.47 = 12.94. 



6.0828/4 = 1.521. Similarly, by comparison with 
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